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ABSTRACT 

We solve for the structure of a hot accretion disc with unsaturated thermal Comptoni- 
sation of soft photons and with advection, generalizing the classical model of Shapiro et 
al. The upper limit on the accretion rate due to advection constrains the luminosity to 
^ 0.15y 3 / 5 a 7 / 5 of the Eddington limit, where y and a are the Compton and viscosity 
parameters, respectively. The characteristic temperature and Thomson optical depth of 
the inner flow at accretion rates within an order of magnitude of that upper limit are 
~ 10 9 K and ~ 1, respectively. The resulting spectra are then in close agreement with 
the X-ray and soft 7-ray spectra from black-hole binaries in the hard state and Seyferts. 
At low accretion rates, bremsstrahlung becomes the dominant radiative process. 

Key words: accretion: accretion discs - black hole physics - galaxies: Seyfert - gamma- 
rays: theory - X-rays: galaxies - X-rays: stars 
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1 INTRODUCTION 

Matter accreting onto a compact object (either in an AGN 
or a binary) is likely to have a net angular momentum, which 
causes formation of a disc. Shakura & Sunyaev (1973) have 
developed a theory of optically-thick disc accretion. Such 
discs emit local radiation close to a blackbody, and thus their 
maximal temperature is at most ~ 10 7 K. Those discs can- 
not thus explain hard X-rays and soft 7-rays observed up to 
~ 1 MeV from both black-hole binaries and Seyfert AGNs. 
Shapiro, Lightman & Eardley (1976, hereafter SLE76) found 
the existence of another branch of disc solution, in which the 
electron temperature, T e , is ~ 10 9 K. In their specific model, 
the dominant radiative process is multiple Compton upscat- 
tering of copious soft photons by the disc electron plasma, 
which forms power-law X-ray spectra extending up to a few 
hundred of keV. 

SLE76 assumed in their analysis that all the dissipated 
energy is locally radiated away. However, advection of hot 
ions is important in a range of parameters of the accre- 
tion flow (Ichimaru 1977). Solutions including both radiative 
cooling and advection have been extensively studied recently 
(e.g., Abramowicz et al. 1995; Narayan & Yi 1995, hereafter 
NY95; Chen et al. 1995; Bjornsson et al. 1996). Some of 
those studies (e.g., Bjornsson et al. 1996) assume Comp- 
tonised bremsstrahlung as the only cooling process, which 
turns out to reproduce well the qualitative properties of the 
solutions. However, the predicted X-ray spectra are then 
very hard, with the photon spectral index of F Ss 1, which 
disagrees with X-ray observations of both black-hole binaries 



and Seyferts, showing typical F ~ 1.7-2 (see, e.g, Zdziarski 
et al. 1997, hereafter Z97, for a review). Other studies (e.g., 
NY95, Narayan 1996) include also Comptonisation of ther- 
mal synchrotron photons in an equipartition magnetic field. 
However, as shown by SLE76, the synchrotron process does 
not supply sufficient flux of soft seed photons to explain the 
X-ray spectrum of Cyg X-l in the hard state. This conclu- 
sion has been extended to GX 339-4 in the hard state and 
Seyferts by Zdziarski et al. (1998, hereafter Z98). Those re- 
sults imply that if Comptonised synchrotron emission were 
the only source of soft photons, the X-ray spectra would be 
harder than those observed. On the other hand, soft-photon 
emission by a cold disc surrounding the hot flow has been 
shown by Esin (1997) either to be negligible or to lead to a 
runaway cooling of the hot flow. 

In spite of these theoretical difficulties, X-ray observa- 
tions show that a supply of soft seed photons sufficient to 
give rise to spectra with F ~ 1.7-2 does exist in luminous 
black hole systems. Also, the shape of a thermal Compton- 
isation spectrum (at energies much above the seed photon 
energy) is a function solely of T e and the Thomson optical 
depth, t. Thus, it is possible to compare a model with X-ray 
data without specifying the soft photon source. Furthermore, 
most of the power in a spectrum with F < 2 is emitted at its 
high-energy end, and then even the model normalization can 
be compared with data without knowledge of the properties 
of soft seed photons. 

Such an approach was adopted by SLE76, who specified 
the local cooling rate in a disc by assuming a constant value 
of the Compton parameter, y, a quantity directly related 
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to the X-ray index F. Observationally, luminous black-hole 
sources typically show an almost constant V with varying 
flux (e.g., Cyg X-l, Gierlihski et al. 1997; GX 339-4, Ueda, 
Ebisawa & Done 1994; IC 4329A, Fiore et al. 1992), which 
justifies that assumption. The origin of soft seed photons 
does not need to be specified in the formalism, but SLE76 
advocated its origin in cold clouds mixed with the hot flow 
(e.g., Kuncic, Celotti & Rees 1997). Thermal synchrotron 
radiation will also provide some soft photons. 

In this work, we generalize the hot disc model of SLE76 
by including advection, by treating the Compton parameter 
as a free paramater (whereas SLE76 assumed y = 1), and 
by extending the solution to low accretion rates. As we show 
below, the model predicts r ~ 1 and T e ~ 10 9 K only weakly 
dependent on the disc radius and the accretion rate. 

Thus, the predictions of the model are in excellent 
agreement with recent studies of broad-band (~ 1-10 3 keV) 
spectra of both black-hole binaries in their hard state and 
Seyferts (e.g., Zdziarski, Johnson & Magdziarz 1996, Z97, 
Z98; Gierlihski et al. 1997; Poutanen et al. 1997; Johnson et 
al. 1997; Magdziarz et al. 1998). Those studies show that the 
intrinsic spectra are well fitted by thermal Comptonisation 
in a plasma with just T e ~ 10 9 K and r ~ 1, with a small dis- 
persion of those parameters between different objects. Note 
that those parameters correspond to y ^ 1. 



2 THE HOT DISC SOLUTION 

In the model of SLE76, gravitational energy is converted into 
ion thermal energy, which is then transferred to electrons by 
Coulomb interactions. The electrons, in turn, radiate their 
energy by Compton upscattering the soft seed photons. 

We assume here a pressure equipartition between the 
magnetic field and gas, P mag = -Pgas, which results in 
P = 2(riikTi + n e kT e ), where Ti, m, T e , n e are the ion and 
electron temperature and number density, respectively. We 
then assume the plasma is two-temperature, i.e., T,: 2> T e 
(which is typically satisfied in an inner part of the disc), 
which allows us to neglect the electron contribution to P. 
The magnetic field here is treated similarly to, e.g., NY95, 
where the viscosity parameter, a, is unrelated to the mag- 
netic field. In another model, Shakura & Sunyaev (1973) 
and Sakimoto & Coroniti (1981) use the magnetic field as a 
source of viscosity, which gives P mag ~ aP sas . Adopting that 
model would slightly modify the equations below, slightly re- 
ducing r for a given accretion rate. We further assume e 
pairs are negligible (as shown to be the case in hot discs 
with advection by Bjornsson et al. 1996 and Kusunose & 
Mineshige 1996), and neglect the radiation pressure. 

We are mostly interested in an inner region of the disc, 
in which a self-similar solution of advection-dominated discs 
found by NY95 breaks down (e.g., Chen, Abramowicz & La- 
sota 1997). Therefore, we consider here local disc solutions 
and parametrize advection using results of global studies, 
following, e.g., Abramowicz et al. (1995) and Bjornsson et al. 
(1996). As those authors, we use the pseudo-Newtonian po- 
tential (Paczyhski & Wiita 1980), $ = -GM/(R-R S ), and 
assume the disc rotation is Keplerian, D. = (GM / R) 1 ^ 2 / (R— 
Rs), where R is the disc radius, Rs = 2GM/<? is the 
Schwarzschild radius, and M is the central mass. We take 
into account the presence of He in the flow. 



With the above assumptions, hydrostatic equilibrium 
and conservation of angular momentum yield (with the 
equation of state of P = 2mkTi), 
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respectively. Here H is the disc scale height, r = Rc 2 /GM, 
q = 1 - 2/r, g = 1 + 4/3rg, J = 1 - (6/r) 1/2 (3g/2) ac- 
counts for the zero-torque boundary condition at r = 6, 
Qi = kTi/rripC 2 , a is the viscosity parameter, r = n e aTH, 
m = M(? /Le, M is the mass accretion rate, the Eddington 
luminosity is Le = 47r/i e GMm p c/<7T, Hi = 4/(1 + 3X) and 
He — 2/(1 + X) are the mean ion and electron molecular 
weight, respectively, X (f» 0.7) is the H mass fraction, and 
<tt is the Thomson cross section. 

The energy is dissipated in the disc at a rate of Q = 
(3/4n)GMMJg/q 2 R 3 per unit area (including both sides 
of the disc). This energy is partly radiated and partly ad- 
vected to smaller radii. We find that a quantity useful to 
describe the local emission of the disc is a local Edding- 
ton ratio, defined as the power per unit lnr in units of Le, 
£ = (dL/dlnR)/L E = 2ttQR 2 /L e . The dissipative £ is, 
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where the function A is needed below. The maximum of £(r) 
occurs at r « 16.9 (compared to r = 27/2 in the Newtonian, 
q = g = 1, J = 1 — 6 1 / 2 /r 1 / 2 , case). A fraction of £ is 
Coulomb-transferred to the electrons and then radiated. The 
rate of Coulomb transfer corresponds to (Spitzer 1956), 
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Here 9 e = kT e /m e c 2 , fa » In A, A w 0.3Q e (H/Ta ) 1/2 , a 
is the Bohr radius, and 9i < e C 1 (which condition is 
typically satisfied in a hot flow). At this limit, the presence 
of He reduces the Coulomb rate with respect to the case of 
pure H by [le/ Hi- This is because the rate is proportional 
to Z 2 /Az, which is unity for both H and He, and ni/n e = 
He/Hi. Here Az is the mass number of a Zth element. When 
the above condition is not satisfied, the rate given for H by 
Stepney & Guilbert (1983) can be generalized to yield an 
fc of, ' 
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where Ki is a modified Bessel function, nz is the density 
of the Zth element, and Q z = Qi/A z . Note that NY95 ne- 
glected the effect of Az in their accounting for the presence 
of He in the Coulomb rate. Therefore, eq. (3.3) in NY95 
and eq. (Al) in Mahadevan (1997) needs to be multiplied 
by 0.8 and 0.8/ie//ii, respectively, in the limit 6; < 6 e C 1; 
otherwise eqs. (|3[)-(U) above should be used. 

SLE76 assumed in their solution that a supply of soft 
seed photons results in the Compton parameter (e.g., Ry- 
bicki & Lightman 1979) of y = 40 e r max(l, r) = 1. This can 
be achieved, e.g., if most of the soft photons are from repro- 
cessing of the hard radiation by cold matter in the vicinity 
of the hot flow (Haardt & Maraschi 1993). A given value of y 
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(which we allow here to differ from 1) approximately corre- 
sponds to a specific value of the energy amplification factor 
of Comptonisation, r\ (at constant soft photon energy, e.g., 
Dermer, Liang & Canfield 1991). Then, the local radiative 
balance can be expressed as, 



^Coul = ri(y)l BO H + ibrcm, 



(•>) 



where £ so { t and Ihtem corresponds to the power in the soft 
photons and Comptonised bremsstrahlung, respectively. For 
now, we interpret ^coui as corresponding to Comptonisation 
of soft photons, which is equivalent to assuming ^coui S> 
4rcm- We also assume Ti ^> T e , which allows us to neglect 
the negative term in eq. (^|). At the end of this Section, we 
discuss the validity of those assumptions and discuss the 
behaviour of the solution when they break down. 

We assume r Ss 1 (which we check it is the case a pos- 
teriori), which corresponds to T e and ^coui of 
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A part of £ is advected to smaller radii (rather than radiated) 
at a rate of Q a d v = (M /2ivR 2 )(P/ Hinim p )£ (Chen & Taam 
1993), where £ is a dimensionless factor found » 1 by global 
hot-disc studies (e.g., Chen et al. 1997). The advected power 
corresponds to, 
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The energy balance of the disc is given by, 

I = ^Goul + ^adv ■ 



(8) 



(9) 



Hereafter, we assume a, f, fc and y to be constant; oth- 
erwise the solutions below of eq. (^) become implicit. The 
assumption of y(r) — constant corresponds to a fixed ratio 
of the soft photon flux to the total local flux (eq. ||), which 
can plausibly result from physical feedback, related to, e.g., 
geometry. The solution has two branches, corresponding to 
dominant cooling (SLE76), i ~ ^Coui tx m, and dominant 
advection, I » £ a dv The first one is, 
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On the other hand, dominant advection corresponds to 
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Note no dependence on a in £coui(r,r). 
The general solution of eq. fl9J) is, 
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where the parameter a = m/r is in the range < a < 
?>ag 2 r x l 2 /2£g (with the lower and upper limits correspond- 
ing to dominant cooling and advection, respectively), and 
the functions A, B, C are defined in equations above. The 
local disc emission corresponds to, 
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Figure 1. Disc solutions at r = 27/2. The left-hand-side curves 
show the optically-thin solution for f = 1, fc = 15. The disc 
emission is due to Comptonisation with y = 1 except for the 
solid curves below the break (at rh ~ 0.002-0.004), for which 
bremsstrahlung dominates. The right-hand-side curves show the 
optically-thick solution for M = IOMq and with a corona above 
the disc dissipating 1/2 of the total power (Svensson & Zdziarski 
1994). Solid curves give the accreted rh whereas the dashed curves 
give m ra d, i.e., the part of rh which is converted into escaping ra- 
diation. Note a good agreement of the m-r dependence presented 
here with those of Abramowicz et al. (1995), Chen et al. (1995) 
and Bjornsson et al. (1996). 



Examples of the solution are shown in Fig. 1 (left-hand- 
side curves). Solid curves give m(r), and dashed curves give 
m ra d(r) = ^coui/^4, i.e., the locally-radiated fraction of m 
(m ra d = rh on the cooling branch). The corresponding tem- 
perature is given by eq. T e ~ 1.5(j//r) x 10 9 K. The max- 
ima of rh, m ra d and r are located in the turning region con- 
necting the cooling and advection branches at a = 3A/4C, 
6A/13C and A/3C, respectively. The maximum m is, 
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Numerically, m n 



e.g., m n 



2.85y 3 / 5 « 7 / 5 J 1 / 5 !? V 7/5 r 3/10 , and, 



- 7.37y i/0 a'' b at r = 27/2. Hereafter, we set 
£ = 1 and fc = 15 (typical for stellar-mass black hole 
sources; for AGNs, fc ~ In A ~ 20) in numerical expres- 
sions. The t corresponding to m max is given by 
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Its maximum of r(m max ) « 1.22y 3/5 a 2/5 is at r = 12.9. 
Advection at m max carries 0.59 of the dissipated energy, and 
the remainder corresponds to ^Coui(?h max ) of, 
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is at 



r = 18.6. The maxima of r and fcoui in the rh-r space 
(Fig. 1) are located in the cooling-dominated region and are 
larger than the values corresponding to m max by the factors 
of (4/3) 4/5 « 1.26 and 2 32 / 5 7 7/5 /13 13/5 « 1.63, respectively. 
We have assumed so far Tj 3> T e , £coui S> 4rem, where 
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^4brem > 1 is a factor including both Comptonisation and 
relativistic corrections and aj is the fine-structure constant 
(Stepney & Guilbert 1983). (Note that the bremsstrahlung 
rate given in eq. [3.5] in NY95 appears incorrect: 1.25 should 
be replaced by /i e .)- The temperature ratio (at r ^ 1) and 
the relative bremsstrahlung rate (at T e <C Tj) are given by, 
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We see that both ratios have very similar dependences on 
the disc parameters. On the advection branch, r oc m, 
and then fbrom/^Coui oc rh~ 2 , which implies bremsstrahlung 
becoming dominant at low rh while still T e <C T. How- 
ever, this happens only at very low ^coui (below the range 
of m ra( j shown in Fig. 1). On the cooling branch, both 
above assumptions break down roughly at the same time. 
Note that the effects on r(m) of the corrections due to 
T e ~ Ti and due to bremsstrahlung are in the opposite 
directions, thus partly canceling each other. Numerically, 
the T e <C T condition (which is marginally more restrictive 
than 4icm <C ^coui) requires rh 3> ya/20 at r = 27/2, and 
r <C 280/(ay) 2//3 at m = 1. Otherwise, the disc is dominated 
by 1-T bremsstrahlung, for which the cooling-dominated so- 
lution is [for P = 2(rii + n e )kT, which modifies eq. ([l])], 
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From eqs. (Q) and (^), we also find that a 1-T disc is char- 
acterized by a specific local y on both branches when r < 1, 
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We see from the above discussion that our main assump- 
tions: T e <C T, ^brcm <C ^Coui, and y(r) — constant, break 
down on the cooling branch roughly simultaneously at low rh 
and/or large r. Thus, in Fig. 1, we show the 2-T, constant-?/, 
solution at large m and the 1-T bremsstrahlung solution at 
small rh. In the latter, we set A^ rBm = 1 since both yir 1 
(eq. |l)) and 6e < 1 (eq. 0) at m < 1. We simply assume 
a sharp transition between the two regimes; in reality, it will 
be smooth and more complex. We note that at radii at which 



bremsstrahlung dominates, m n 



-1/2 



(Abramowicz et 



al. 1995), whereas rh max oc r 3//1 ° in the constant-y, 2-T, re- 
gion (eq. Q). 

In contrast to the above results, Liang & Thompson 
(1979) and Wandel & Liang (1991) extended the 2-T so- 
lution of SLE76 (eq. Jl(]]) to low values of rh obtaining a 
cooling-dominated, 1-T solution with y = constant, rh oc 
t~ , t > 1. That extension neglects bremsstrahlung; how- 
ever, we find the bremsstrahlung power on that solution ex- 
ceeds the dissipated power, with the ratio of the two oc r 3 . 
Thus, that solution branch appears not to exist. Instead, the 
behaviour of the SLE76 solution at low rh is given by a tran- 
sition to the 1-T bremsstrahlung model, as shown above. 

Note that y is necessarily variable in 1-T flows with 
r < 1 (eq. |Fl]]) as well as it becomes large at small r, e.g., 
y X T ~ 7rh/a at r — 27/2. This represents a strong argument 
against their presence in inner regions of luminous black-hole 
systems, with constant y < 1 observed (Section 1). 



3 DISCUSSION 

We have solved the local structure of a hot accretion disc 
with the main cooling process being thermal Comptonisa- 
tion of soft photons. Some feedback processes are assumed 
to lead to y — constant in a range of radii, which corre- 
sponds to an approximately constant X-ray spectral index, 
r, of emission from each radius. The total spectrum can be 
then obtained by radial integration. The accretion rate of the 
hot flow can be constant throughout the flow, or, plausibly, 
it can equal the local maximum rate (eq. jla|). The latter 
case corresponds to a a two-phase flow, in which a cold flow 
(e.g., in the form of clouds) carries the difference between 
the total M and that corresponding to the maximum M in 
the hot flow. The latter case maximizes the luminosity, 



, n T r 3/5 7/5 r 

1 0.15?/ a L EM , 



(22) 



which is obtained by radial integration of eq. (|lj]) (0.5Lhot 
comes then from r 40 and the upper limit of integration is 
not essential as long as it is ^> 40). That case also yields the 
plasma parameters only weakly changing with r, in particu- 
lar r ~ 1 over the hot flow (eq. |H|). Interestingly, the limit 
of eq. ( p^ ) appears to agree with the maximum luminosities 
observed from black-hole X-ray novae and persistent X-ray 
binaries in the hard state (Tanaka & Lewin 1995; Tanaka 
& Shibazaki 1996; Phlips et al. 1996; Harmon et al. 1994) 
and AGNs (Peterson 1997) provided a is large (cf. Narayan 
1996; Z98). 

The optical depth of the hot flow close to M max is r ~ 1 
(eq. (Tell), and T e ~ 10 9 K (eq. §), which agree well with 
the observed X-ray spectra (see Section 1). The values of r 
and T e are basically independent of M (except for a slight 
dependence via the Coulomb logarithm), which can explain 
the close similarity of the broad-band spectra of black-hole 
binaries in the hard state and Seyferts (e.g., Z97, Z98). Also, 
the solutions very weakly depend on the advection param- 
eter, £ [e.g., r(m max ) oc £ -1//s ], which accurate value is ob- 
tainable only from global studies. In fact, a more important 
effect is geometry of Comptonisation models, which should 
be specified in comparing fit results with hot disc models. 
E.g., fits in Z97 are for a hemisphere geometry (Poutanen 
& Svensson 1996), in which case the fitted radial r is higher 
than the half-thickness r obtained from fitting a slab model. 
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For example, Z98 obtains for GX 339-4 r ~ 2 and 0.8 in the 
hemisphere and slab geometry, respectively. The latter ge- 
ometry is more appropriate for hot discs radiating close to 
M max , with relatively important cooling. 

Black hole systems (in the hard state) show weak spec- 
tral variability of the spectral shape with changing flux (see 
Section 1). We can explain this by the weak dependence of r 
(determining the spectral shape) on ^ooul- This dependence 
is t oc £q^ u1 on the cooling-dominated branch, and r tx ^coul 
on the advection-dominated branch (eqs. @-@). Which 
branch corresponds to existing accretion discs can thus be 
determined by studying the X-ray spectral variability. 

The cooling-dominated branch is viscously stable, but 
unstable thermally (Pringle 1976). We note, however, that 
this thermal instability is different and much weaker from 
the classical bremsstrahlung-cooling instability of Pringle, 
Rees & Pacholczyk (1973), in which an increase of tempera- 
ture results in a decrease of the cooling rate per unit area. In 
the model of SLE76, an increase of Ti yields an increase of 
the cooling rate, but the corresponding increase of the local 
heating is faster. The flow becomes stable when advection 
is significant (e.g., Wu 1997). 

Finally, we compare the hot solution with the optically- 
thick (cold) disc solution of Shakura & Sunyaev (1973) and 
Abramowicz et al. (1988). The right-hand-side curves in Fig. 
1 show the latter solution, but with 0.5 of the dissipation 
occuring in a corona (Svensson & Zdziarski 1994). The solu- 
tion above m ~ 1 (weakly dependent on a) has Am /At < 
and is both viscously and thermally unstable (Pringle 1981). 
(The flow becomes again stable due to advection of radiation 
when L ~ Le, Abramowicz et al. 1988, with Le correspond- 
ing to rh = 16). Thus, the maximum stable m ~ 1 of the 
cold solution in Fig. 1 corresponds to ~ 5% of Le- Notably, 
this coincides with the luminosity of Cyg X-l in the soft 
state (Gierlihski et al. 1998), which supports identification 
of that state with the lower turning point of the cold solu- 
tion. (The presence of a corona, assumed above, is inferred 
from observations of power law tails extending beyond the 
disc blackbody in the soft state of black-hole binaries, Grove 
et al. 1998; Gierlihski et al. 1998.) We see in Fig. 1 that the 
fhmax of the hot solution can be, depending on its a, either 
larger or smaller than the maximum stable m of the cold 
flow. Thus, the soft-to-hard state transition in X-ray bina- 
ries (e.g., Tanaka & Lewin 1995) can in principle correspond 
to either an increase or decrease of M. 



ACKNOWLEDGEMENTS 

This research has been supported in part by the Polish 
KBN grants 2P03D01008 and 2P03D00614. I thank Marek 
Gierlihski for valuable comments on this work and help with 
testing the solutions, and Bozena Czerny for discussions and 
supplying her optically-thick disc code. 

REFERENCES 

Abramowicz M. A., Chen X., Kato S., Lasota J. -P., Regev O., 

1995, ApJ, 438, L37 
Abramowicz M. A., Czerny B., Lasota J. -P., Szuszkiewicz E., 

1988, ApJ, 332, 646 



Bjornsson C, Abramowicz M. A., Chen X., Lasota J. -P., 1996, 
ApJ, 467, 99 

Chen X., Abramowicz M. A., Lasota J. P., Narayan R., Yi I., 

1995, ApJ, 443, L61 
Chen X., Abramowicz M. A., Lasota J. P., 1997, ApJ, 476, 61 
Chen X., Taam R. E., 1993, ApJ, 412, 254 
Dermer C. D., Liang E. P, Canfield E., 1991, ApJ, 369, 410 
Esin A. A., 1997, ApJ, 482, 400 

Fiore F., Perola G. C, Matsuoka M., Yamauchi M., Piro L., 1992, 
A&A, 262, 37 

Gierlihski M., Zdziarski A. A, Coppi P. S., Poutanen J., Ebisawa 
K., Johnson W. NT., 1998, Nucl. Phys. B, in press 

Gierlihski M., Zdziarski A. A, Done C, Johnson W. N., Ebisawa 
K., Ucda Y., Haardt F., Phlips B. F., 1997, MNRAS, 288, 958 

Grove J. E., Johnson W. N., Kroeger R. A., McNaron-Brown K., 
Skibo J. G, 1998, ApJ, 499, in press 

Haardt F., Maraschi L., 1993, ApJ, 413, 507 

Harmon B. A., et al., 1994, ApJ, 425, L17 

Ichimaru S., 1977, ApJ, 214, 840 

Johnson W. N., McNaron-Brown K., Kurfess J. D., Zdziarski A. 
A., Magdziarz P., Gehrels N., 1997, ApJ, 482, 173 

Kuncic Z., Celotti A., Rees, M. J., 1997, MNRAS, 284, 717 

Kusunose M., Mineshige S., 1996, ApJ, 468, 330 

Liang E. P. T., Thompson K. A., 1979, MNRAS, 189, 421 

Magdziarz P., Blaes O. M., Zdziarski A. A., Johnson W. N., Smith 
D. A., 1998, MNRAS, in press 

Mahadevan R., 1997, ApJ, 477, 585 

Narayan R., 1996, ApJ, 462, 136 

Narayan R., Yi I., 1995, ApJ, 452, 710 (NY95) 

Paczyhski B., Wiita P. J., 1980, A&A, 88, 23 

Peterson B. M., 1997, Active Galactic Nuclei, Cambridge Univer- 
sity Press, Cambridge 

Phlips B. F., et al., 1996, ApJ, 465, 907 

Poutanen J., Svensson R., 1996, ApJ, 470, 249 

Poutanen J., Svensson R., Stern B., 1997, in Winkler C, Cour- 
voisier T., Durouchoux P., eds., ESA SP-382, The Transparent 
Universe, 401 

Pringle J. E., 1976, MNRAS, 177, 65 

Pringle J. E., 1981, ARA&A, 19, 137 

Pringle J. E., Rees M. J., Pacholczyk A. G., 1973, A&A, 29, 179 
Rybicki G. R., Lightman A. P., 1979, Radiative Processes in As- 
trophysics, Wiley-Interscience, New York 
Sakimoto P. J., Coroniti F. V., 1981, ApJ, 247, 19 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Shapiro S. L., Lightman A. P., Eardley D. M., 1976, ApJ, 204, 
187 (SLE76) 

Spitzer, L., 1956, Physics of Fully Ionized Gases, Interscience, 
New York 

Stepney S., Guilbert P. W., 1983, 204, 1269 
Svensson R., Zdziarski A. A., 1994, ApJ, 436, 599 
Tanaka Y., Lewin W. H. G., 1995, in: Lewin W. H. G., van 
Paradijs J., van den Heuvel E. P. J., eds., X-Ray Binaries, 
Cambridge University Press, Cambridge, 126 
Tanaka Y., Shibazaki N., 1996, ARA&A, 34, 607 
Ueda Y., Ebisawa K., Done C, 1994, PASJ, 46, 107 
Wandel A., Liang E. P., 1991, ApJ, 380, 84 
Wu X.-B., 1997, MNRAS, 292, 113 

Zdziarski A. A., Johnson W. N., Magdziarz P., 1996, MNRAS, 
283, 193 

Zdziarski A. A., Johnson W. N., Poutanen J., Magdziarz P., 
Gierlihski M., 1997, in Winkler C, Courvoisier T., Durou- 
choux P., eds., ESA SP-382, The Transparent Universe, 373 
(Z97) 

Zdziarski A. A., Poutanen J., Mikolajewska J., Gierlihski M., Ebi- 
sawa K., Johnson W. N., 1998, MNRAS, submitted (Z98) 



